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Abstract 

We analyse the possible presence of New Physics (NP) in the 
Flavour Sector and evaluate its potential for solving the tension be- 
tween the experimental values of Aj/^k^ and Br(S+ t^i^t) with 
respect to the Standard Model (SM) expectations. Updated model 
independent analyses, where NP contributions are allowed in S^-S^ 
and Bg-Bg transitions, suggest the need of New Physics in the bd 
sector. A detailed analysis of recent Flavour data is then presented 
in the framework of a simple extension of the SM, where a Q = 2/3 
vector-like isosinglet quark is added to the spectrum of the SM. Spe- 
cial emphasis is given to the implications of this model for correlations 
among various measurable quantities. We include constraints from all 
the relevant quark flavour sectors and give precise predictions for se- 
lected rare processes. We flnd important deviations from the SM in 
observables in the bd sector like the semileptonic asymmetry ^f^, 
B^ fi^ li~ and — ^f^. Other potential places where NP can 
show up include 7, n^vv^ t ^ Zq and ji^ ji~ 

among others. The experimental data favours in this model the exis- 
tence of an up vector-like quark with a mass below 600(1000) GeV at 
1(2) a. 

^ fbotella@uv.es, ^gbranco@ist. utl.pt, ^negomi@uv.es 



1 Introduction 



The flavour puzzle remains one of the fundamental questions in particle 
physics |1|. In the Standard Model (SM), the flavour structure of Yukawa 
couplings is not constrained by gauge invariance, which in turn implies that 
quark masses and their mixing are free parameters of the theory. Apart from 
this important shortcoming, the SM has been very successful in describing 
flavour mixing and CP violation through a 3 x 3 unitary Cabibbo-Kobayashi- 
Maskawa (CKM) |2^^ matrix characterized by four independent parameters. 
At present, the SM and its built-in CKM mechanism for mixing and CP vi- 
olation are in good agreement with most of the experimental data. This 
is an impressive achievement of the SM, given the large amount of data. 
The proliferation of free parameters in the SM is one of the motivations for 
considering New Physics (NP) which could shed some light on the flavour 
puzzle. On the experimental side, there is also motivation to consider NP 
with flavour implications; in particular, there are hints of potential devia- 
tions from SM predictions, such as: 

(i) A tension between the experimental value of ^j/^k^ j4|[5] (the time 
dependent CP asymmetry in J/ipKs) and the value implied by flavour 
flts. In the context of the SM, ^j/^^^ = sin 2/3 where /3 = arg [— V^^V^^KdKt]- 
The experimental value of ^ j/^x^ is around 2a lower than the value of sin 2/3 
extracted indirectly from other experimental input. 

(ii) The experimental value of the branching ratio Br(S+ t^i^t) is 
around 2.5a larger than the value expected in the SM (6]-[9]. 

In this paper, we examine the question whether the present data may 
already give some hints of New Physics (NP) and also discuss the potential 
for data expected in the future to give conflrmation and /or provide more 
restrictive bounds on NP. We start our analysis by presenting a brief update 
of a model independent analysis (MIA) looking for the presence of NP [lO,ll 



in view of the present experimental data. For deflniteness, in the MIA, we 
shall make the following assumptions: 

(i) We assume that the extraction of parameters from SM tree level dom- 
inated processes is not affected by the possible presence of NP, like in: 
Kdl \VJ, \VJ, \VJ, Ksl \VJ and 7 = arg [-^.K^KT^K*.]- 

(ii) We allow for NP to give signiflcant contributions to B^-B^ and B^-B^ 
transitions. 

In the MIA we shall assume that the CKM matrix 1/ is a 3 x 3 unitary matrix. 
Our MIA applies to a large class of theories beyond the SM, including super- 
symmetric extensions of the SM and some generic two Higgs doublet models. 
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However, it should be emphasized that 3x3 unitarity is an assumption which 
can be violated in reasonable extensions of the SM. 

The main goal of this paper is to make a study of the impact of NP 
on the flavour sector, in the context of a simple extension of the SM where 
there are naturally small violations of 3 x 3 unitarity in V. For deflniteness, 
we consider an extension of the SM where an isosinglet (5 = 2/3 vector-like 
quark is added to the SM. The paper is organised as follow. First we present 
the updated MIA. In section |3] we deflne the model with a singlet up vector- 
like quark (sUVL), then we successively present the flavour and CP results 
in the bd^ bs and ds sectors. We devote special attention to correlations 
among different observables. In section [7| we discuss rare up decays and the 
mass of the new up quark. Then we present our conclusions and devote two 
appendices to explain additional details. 



2 Model Independent Analysis 

With the standard deflnitions of the rephasing invariant quantities [12 



a 


= arg 


[-Vt,v,,v,iv:,] 


/3 


= arg 
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= arg 


\—V V V* V* 

I ^ ud^ cb^ ub^ cd 




= arg 





J y 



(1) 

a + ;5 + 7 = 7rby deflnition, independently of whether the CKM matrix V 
is 3 X 3 unitary or not. It is well known that in the context of the SM, one 
has, to an excellent approximation: 



"""Z*^^ sin {AMs, t) [r (5° ^ J/i,Ks) + T (5° ^ J/i,Ks)] ^ ' 

(2) 

In a similar way, in the SM, the time-dependent asymmetry in {7T7t)j=2 
vs. {^^)i=2'j ^(7r7r)/^2 5 whcrc {7111)1=2 denotes the strong isospin state 

of the pair of pions, is 

^(7r7r)j^2 = sin(2a) . (3) 
As 2a = — (;5 + 7) mod [27r], A^^^^^)^^^ can be viewed as a measurement of 



7; direct measurements of 7 are achieved, anyway, in other channels 13 -20 
In the SM, Aj/^x^, ^(7r7r)j^2 measurements and moduli of \V-j\ are related 
through the unitarity relation 



^udY^ub 



sin(7 + /5) = sin(/3), (4) 
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which constitutes a nice consistency check of the SM flavour picture. 

When new physics is considered, a minimal and reasonable assumption is 
to allow for new contributions to loop-controlled processes but not for tree 
level observables, while preserving 3x3 unitarity of the mixing matrix. In 
particular, potentially modifled mixings of B mesons (both and Bs) can 
be written in the following form 
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SM 



(5) 



where 



J SM 



stands for the SM contribution and {^^,0^} parameterise 

NP-induced deviations from SM expectations. The CP asymmetries in eq.([2]) 
and eq.([3]) are automatically modifled to 

Aj/^ifKs = sin(2(/3-(/)rf)) = sin(2^) , ^(7r7r)j^2 = sin(2(a + (/)rf)) = sin(2a) . 

(6) 

It is clear that 7, or equivalent ly a + /3, can be extracted from experiment in 
a model independent way, as 



TT — 7 = (arcsin Aj/^Xs + arcsin74(^^)^^2) • 



(7) 



Then, using eq.(4) with Ry, = 



tan ( 



, one can obtain 



Ru sin(7 + /3) — sin /5 
cos ^ — Ru cos(7 + ^) 



21 



(8) 



In the absence of NP, f3 = f3 and (^d = 0. Using eq.([8|) and taking into 
account experimental data, one obtains 



tan ( 



0.11 ±0.03 



A complete analysis yields the result summarized in flgure [TJ 

The conclusion is clear: there is a tension between i?^, ^ and a, and NP 
in B^-B^ may solve or relax it. Although dominated by the SM picture, 
data appears to be pointing to a small but signiflcant presence of NP in the 
bd sector. 

One can extend the analysis to the bs sector | TH[22] - [25] by considering 
the time-dependent asymmetry (contrary to B^ J/^ Ks-> angular analysis 
to separate different CP components in the flnal state is required) in 5^ 
vs. 5° Aji^^, which in the SM is 



= sin(2/35) , 



(9) 



4 




expected to be small, since 3x3 unitarity, in the CKM framework, forces 
f3s = 0.187 ± 0.006. Allowing for the modification in eq.([5]) gives 

^j/^^ = sin(2(/3, + (/),)). (10) 

It is, therefore, the LHCb measurement of Aj/\^^ - and related channels - 



that dominates the extraction of (f)s [26,27 . There is ample room for NP in 
the bs sector - (j)s can be two or three times larger than f3s However, the 
current situation, dominated by experimental uncertainty while being close 
to SM expectations does not require the presence of NP. 



3 Up vector-like singlet quark model 

For definiteness, we consider an extension of the SM where one isosinglet 
vector-like quark T with charge Q = 2/3 is added to the spectrum of the 



SM 28-54 . The 3x3 mixing matrix connecting standard quarks is no 
longer unitary, but a submatrix of a 4 x 4 unitary matrix U. Without loss of 
generality, one can choose to work in the weak basis (WB) where the down 
quark mass matrix is diagonal and real. In this basis, [/ is a 4 x 4 unitary 
matrix which enters the bidiagonalization of the 4x4 mass matrix A4 of 
(5 = 2/3 quarks: 

U^MM^U = diag. {ml, ml m\, mT^) . (11) 
The charged and neutral current interactions have the form 

= -^Ul 1^ VdLW, + h.c. , (12) 

= -TT^ir ["^ ^'^ (^^^) -dLrdL-2 sm2 JL] , (13) 

Z COS u\Y 
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where d = {d, s,b), u = {u, c, t, T) and \^ is a 4 x 3 submatrix of the matrix 
U: 



U 
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Vus 






\ 




Vcs 










Vts 




Uu 




\ Vrd 


Vts 






J 



(14) 



It is clear that the submatrix V(3x3), i.e. the upper left 3x3 block within 
[/, is not a unitary matrix, since ^^3^3) ^(3x3) 7^ 1(3x3)- However, these devi- 
ations of unitarity of the "would-be standard" mixing matrix are naturally 
suppressed by the ratio m^/m|., where m denotes generically the standard 



quark masses 28-31]. These deviations from unitarity lead to flavour chang- 
ing neutral currents (FCNC) which are present just in the up sector and 
controlled by 



Sij - u. 



ui 



(15) 



These FCNC are thus suppressed by the ratio /m\. This natural suppres- 
sion of FCNC is crucial in order to make the model plausible, since FCNC 
mediate dangerous tree-level processes leading for example to D^ D^ mixing 
though [/^4t/*4. 



3.1 Comparison with the 4-generations Standard model 

At this stage, it is instructive to compare the sUVL model with the four gen- 
erations Standard Model (SM4). The differences between these two minimal 
extensions of the SM can be summarized as follows: 



Both the sUVL model and the SM4 involve deviations of 3 x 3 uni- 
tarity in V . In the sUVL model deviations of unitarity are naturally 
suppressed by the ratio /m\ and lead to naturally small Z-mediated 
FCNC in the up sector. In the SM4 there are no tree level FCNC due 



to an exact GIM mechanism 55 



2. The dominant Higgs boson production mechanism at the LHC is the 
gluon-gluon fusion process which is dominated by a heavy quark loop 



56,57|. With the two new heavy sequential quarks introduced in the 
SM4, there is the potential of having an increase in the production cross 
section by a factor of (1 + 2)^ = 9. Even if smaller, this enhancement 
has very important consequences [58]-j61j . In the sUVL model there is 
no enhancement, since the quark mass terms of the vector-like quark 
do not arise from electroweak symmetry breaking. 
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4 The bd sector 



As a starter, let us analyse how the sUVL model is able to mitigate the 
tensions present in the bd sector of the SM. Figure [2] displays the usual 
unitarity triangle (SM) together with a realistic quadrangle corresponding 
to the sUVL model. It is important to stress that, in order to reduce the 
tension, we need an increase of to enlarge Br(5+ t^i^t)- But, in the 
SM, this comes hand in hand with an enlarged (5^ thus generating a conflict, 
a tension, with the experimental value [62^|63j . This tension can be solved by 
the necessary NP in the mixing, which, in this model, is flxed by the enlarged 
mixing matrix and the mass of the heavy quark. Focusing on the hd sector, 
this means that it is flxed by the shape of the hd unitarity quadrangle. So, in 
practice, one has more freedom by increasing \V^^\^ increasing /5 and adding 
the fourth side IV^^V^^I, necessary to match the potential modiflcations in 
\VJ and /3. 




^cdYch ^cdYch 

Figure 2: Unitarity quadrangle vs. unitarity triangle. 



A particular case of how the NP in the mixing transforms (3 into 
corresponding to the previous quadrangle, is shown in the flgure [3) Notice 
how the SM term ^ the one to the left, oriented with respect to the horizontal 
(KdKl)^ corresponding to the box diagram with internal t quarks, receives 
additional contributions from the box diagrams with both T and t and only 
with T, respectively. As one can see, starting from a large 2(5 one ends up 
with a smaller 2^. The diflFerent S'o functions are the well-known Inami-Lim 
functions (64). In order to show that the sUVL model certainly solves the 
tensions in the hd sector we flrst present together the results of the SM- 
CKM and sUVL flts in the (Br (5+ ^ r+z/^), ^j/^^J plane in flgure Q The 
experimental input used, together with the methodology and further details, 
are explained in the appendices. Let us just clarify here that, in flgure [3] (as 
in many other flgures along the present paper), are represented three (blue) 
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regions corresponding to 68%, 95% and 99% CL (darker to lighter). The 
eUipses represent the same CL regions for experimental data. 




Figure 3: Pictorial representation of the different contributions to M]^2'^ 
the sUVL model. 



Figure 4(a) clearly shows the tensions among these observables in the SM 
case. In figure 4(b)| it is easy to appreciate how the sUVL model solves the 
tensions in the bd sector. Of course, in order to obtain in this model the 




12. 



10. 



0.64 0.66 0.68 0.7 0.72 0.74 0.76 
(a) SM case. 



22.5- 
20.0- 
17.5- 
15.0 
12.5- 
10.0 
7.5: 




0.64 0.66 0.68 0.7 0.72 0.74 
(b) sUVL model case. 



Figure 4: Br(S+ ^ r+z/^) x 10^ vs. Aj/^Ks^ 68%, 95% and 99% CL regions 
(darker to lighter); the ellipses show the 68%, 95% and 99% CL regions 
corresponding to the measurement. 
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experimental value of Br (5+ r^v^)^ the only possibility is through an 
increase of IV^^I, as can be seen in figure 



changing ^j/^k^, as shown in figure 5(b) 



^ure 


5(a)| 


5(b) 






3.5 4.0 4.5 5.0 
(a) iKbl X 10^. 




0.60 



0.65 0.70 

(b) Aj/^xs- 



0.75 



Figure 5: Ax^ profiles of \Vuh\ and Aj/x^jKs sUVL model. 



A first consequence is that this can be done by replacing the bd triangle 
by a bd quadrangle whose fourth side is given by l^^^^^l- 



0.04 




01 ^ ^ ^ ^ ^ 

0.05 0.1 0.15 0.2 0.25 

Figure 6: \VTd\ vs.|l^5|, 68%, 95% and 99% CL regions (darker to lighter) in 
the sUVL model. 



Figure [6] shows that a typical value of the fourth side of the quadrangle 
is IV^^V^^I ^ A"^. This fourth side allows to enlarge and to modify /3 
in such a way that the final value of ^ yields Aj/\^Xs perfect agreement 
with the experimental value. Following [iT], it is convenient to write down 
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the semileptonic asymmetry for in terms of physical observables. 



A 



SL 



AM, 



[{b + c-a)\V^,Vj\m{2a) + 



+ (2c - a) \V,aVuiycdVc,\ sm (2^ + 7) - c \V,,Vj Aj,^Ks] ■ (16) 

It has to be stressed here that we have included the NP that this model 
introduces in M^^ ^ but the model does not introduce any new contribution 
in ' Of course, one must be careful in not using any expression where 



The first thing to note in eq.(16) is that 



3x3 unitarity has been used [65 -69 
the quadratic term in \V^^\ has a coefficient sin {2a) = 0.00 ± 0.15 so, at the 
leading order, A^j^ is linear in IV^^I- It was also noted in flT^ that there is a 
strong cancellation among the second and third terms so one can guess that 
a 10-20% enhancement in can be easily translated into an enhancement 
factor of five or more in A^j^. This can be easily in seen in figure [7| 




7.5 10.0 12.5 15.0 17.5 20.0 22.5 



Figure 7: A^sl x 10^ vs. Br(5+ ^ r+z/^) x 10^ 68%, 95% and 99% CL regions 
(darker to lighter); the ellipses show the 68%, 95% and 99% CL regions in 
the SM fit. 



This opens the door to an important contribution to the dimuon asym- 
metry observed by the DO collaboration [70 j. We will come back to this point 
in the sequel. Nevertheless, it has to be pointed out that this is a common 
feature of any model where one enlarges \V^i^\ and introduces NP in the B^- 

mixing in order to compensate the potential mismatch of ^j/^k^ coming 
from a /3 different from ^. 

To have a clear picture of the most important changes in the unitarity 
triangle - now a quadrangle - we plot /3 in the SM and in the sUVL model 
fits in figures 8(a) and |8(b)| respectively. 
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(a) SM case. (b) sUVL case. 

Figure 8: Ax^ profile of the phase /3. 



Another important consequence of having a relatively large value of {Vj^^Vj^ 
is an important impact in b ^ d transitions like in the process 
We plot, in figure [9| the likelihood profiles of Bi^Bd both in the 

SM and in the sUVL model. 




0.8 1.0 1.2 1.4 1.6 1.8 5 10 15 20 

(a) SM case. (b) sUVL case. 

Figure 9: Ax^ profile of Br(Srf M^M") x 10^^. 



It is clear that a large enhancement - a factor of five to ten - still is 
possible in this process. It is remarkable that the last LHCb upper bound 
is the most important constraint on this process 71-75|. In the future 
7r^/x^/i~, for which LHCb recently announced a preliminary re- 
sult, will also have a role to play. From the correlation among Bi^B^ 

and Br(5+ ^^^r), shown in figure 10(a), it is clear that keeping 
Br (5+ t^i^t) around its actual experimental central value implies that 
Br(5(i has to be seen at a rate larger than the SM value. 
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A similar correlation can be observed among large negative A^j^ values 
and having BT{Bd larger than the SM value. This can be seen in 

figure [TO(E)| 




7.5 10.0 12.5 15.0 17.5 20.0 22.5 2.5 5.0 7.5 10.0 12.5 15.0 

(a) BT{Bd X 10^° vs. (b) x 10^ vs. BT{Bd x 

Br(5+ ^ T+Ur) X 10^ 10^^ 



Figure 10: Correlations in the sUVL model, 68%, 95% and 99% CL regions 
(darker to lighter) are shown. 
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5 The bs sector 



Last summer, LHCb announced, at the Lepton-Photon Conference, its first 
measurement of the CP asymmetry to final CP eigenstates in the 

decay J/ip^ [26,27 . The associated value of the CP violating phase, 

2(3 s in the SM, was restricted to be small: this was qualified by some speaker 
as the Lepton-Photon drama. This designation can be understood looking 



at the correlation plot |ll(a) , where we can see the 68%, 95% and 99% CL 
regions in the plane {Aj/^,^^ A^sl}-> with A^^j^ the semileptonic asymmetry in 
decays. 



7.5 



5.0 



2.5 



0.0 



-2.5- 



-5.0 



-7.5 



-0.1 0.0 0.1 0.2 

(a) A%^x 10"^ vs. Aji^^. 



0.10 



0.05 



0.00 



-0.05 



-0.10 




-0.1 0.0 0.1 

(b) 13s vs. Aji^^. 



0.2 



Figure 11: Correlations in the sUVL model involving Ajjy^^] 68%, 95% and 
99% CL regions (darker to lighter) are shown. 



If the large DO dimuon asymmetry had to be explained mainly 

through a large coming from NP in the hs sector, this would imply 
also a large ^j/^<|), in confiict with the LHCb result. The origin of the name 
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As one can see in figure 



76 


77 


were 


figure 


11(a), 



compatible with large values of Aj/y^^ 
once the LHCb results on Aj/y^^ are included, A^^j^ is constrained to be in 
the 10~^ range. 

In spite of the fact that there are not very large weak phases in the 
B^^B^ mixing, from the previous figure we must stress that values of A^^^ 
and Aj/y^^ much larger than the SM value are allowed. These values are 
completely correlated. Another important correlation pointed out in |52j is 
among Aj/y^^ and 13 g- Note that in the SM this correlation - for the chosen CP 
final eigenstate - is positive, but in the sUVL model it is negative. This was 
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proven in reference [52] and therefore this correlation provides a consistency 
check of the entire analysis. It is shown in figure ll(b)[ In figure [12] we 
see that after the impressive bounds on the branching ratio Bt{Bs 
presented by LHCb |72|, CMS [74j and ATLAS [75j, there is stiU room for 
values of Aj/x^,^ of order 0.1 . 



0.2 



0.1 



0.0 



-0.1 




4. 



Figure 12: Aj/^^ vs. Br(5, x 10^ in the sUVL model (68%, 95% 

and 99% CL regions, darker to lighter). 



It is apparent that smaller values - below the SM value - of Bt{Bs 
are more easily related to important deviations of Aj/^<|) from its SM 



value, especially towards 0.1 . In fig. 13 we can see two important results. 



First, the model cannot reproduce the value of reported by the DO 
collaboration, in spite of the fact that the model violates 3x3 unitarity, a new 
ingredient generally not included in the majority of the model independent 
analysis. The maximum value allowed in this model for ^4^^ is around — 2 x 
10~^, a factor of seven larger than the SM value but still at 3a from the DO 
reported value. Second, we also learn that the leading contribution in these 
models to the dimuon charge asymmetry comes from Af^, as it should after 
the Lepton-Photon drama which implies a very strong indirect constraint on 



"^SL- completeness we also give, in fig. 13(b), the correlation among 
and the quantity which is expected that LHCb could deliver soon: A^^j^ — A^j^. 

The dominance of A^j^ versus A^^j^ is again clear. In order to clarify 
the impact of future improvements on the measurement of 7, we show the 



correlation between Aj/y^f^ and 7 in figure 14, It is clear that the largest 
departures of Aj/\^^ from the SM value can be easily accommodated in the 
region with larger 7. 
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-2 



-3 



-4 
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-4 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 
(a) A%j^ X 10^ vs. X 10^. 



12 3 4 

(b) A^^xlO^ vs. (A|^-A^^)xl03. 



Figure 13: Correlations with in the sUVL model, 68%, 95% and 99% 
CL regions (darker to lighter). 



Interesting correlations in this model exist among ^i{Bs 1^^ and 
Br (5 -^sT) or Br (5 Xsji^ It is important to stress that low values 
of ^r{Bs 1^'^ are anticorrelated with large values of Br (5 -^sT) and 



Br(5 XsjJi^ jJi ) as can be seen in fig. 15 



0.2 



0.1 



0.0 




-0.1 

1. 1.1 1.2 1.3 1.4 

Figure 14: Aj/s^j^ vs. 7 in the sUVL model, 68%, 95% and 99% CL regions 
(darker to lighter). 
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6. 



4. 



2. - 
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3.3 3.4 3.5 

(a) Bt{Bs 10^ vs. Br(5 




1.5 2. 2.5 3. 

(b) Br(55 ^ l^^fi~)x 10^ vs. Br(5 ^ 



Figure 15: Correlations with Bi^Bg ji^ ji ) in the sUVL model, 68%, 95% 
and 99% CL regions (darker to lighter). 



6 The ds sector 

Of course we have included the most restrictive kaon constraints in the anal- 
ysis, both CP violating and CP conserving |78|. As we have seen in previous 
sections, to eliminate the tensions in the hd sector we need non negligible 
Vj.^ and Vrj.^. It seems that we can have significant NP contributions to 
Br(X+ TT^jyiy) provided we have important It is worthwhile 

to look to the correlation among Br(X+ n^uu) and IV^^I, but also with 
^j/^<^, that can be sizeable provided one has a significant bs quadrangle. 

Both figures 16(a) and |16(b)] confirm our expectations. Bt{K^ n^uiy) 
clearly has a branch that grows with and we can see that positive 

values of Aj/x^f^ - larger than the SM value - can accommodate values of the 
Bt{K^ n^uu) enhanced with respect to the SM values by factors of two 
or more. It is important to point out that in this model the Grossman-Nir 
bound [79] in the plane Bt{K^ Tr^iyu) vs. Br(i^'i, tt^uu) cannot be 
saturated even if the process tt^z/z/ can be enhanced by factors up to 



seven as it is shown in figure 17 
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0.01 0.02 0.03 0.04 -0.1 0.0 0.1 0.2 

(a) Br(i^+ ^ TT^uu) x 10^^ vs. \Vts\^ (b) Br(i^+ ^ tt^uu) x 10^° vs. Aj/^<|,. 



Figure 16: Correlations with Br(X+ n^uu) x lO-*^^ in the sUVL model, 
68%, 95% and 99% CL regions (darker to lighter). 



3.0 




0.5 1.0 1.5 2.0 2.5 3.0 



Figure 17: Bt{Kl ^ ttVp) x 10^° vs. Br(X+ ^ n^uu) x 10^^ in the sUVL 
model, 68%, 95% and 99% CL regions (darker to lighter). 
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7 The up sector 



7.1 Rare top decays 



A very distinctive signal of this class of models are the rare top decays. In 
this case, due to the presence of flavour changing neutral currents at tree 
level in the up sector, we have mainly t ^ cZ and t uZ controlled by 
\UcAUtA\ and \Uu4Ut4\ respectively. Its natural order of magnitude is at the 
10~^ level arriving in some cases to 10~^ as can be seen in the correlation 



plot in flgure 18 
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Figure 18: Br(t cZ) x 10^ vs. Br(t uZ) x 10^ in the sUVL model, 
95% and 99% CL regions (darker to lighter). 



We also display the individual likelihood proflles to see clearly their nat- 
ural order of magnitude. With the values shown in flgure 19, although very 
large compared to the SM values, it is going to be diflicult to discover this 
class of NP using these rare top decays. 



7.2 Charm 

With the presence of FCNC at tree level in the up sector, the D^-D^ mixing 
can be explained naturally in the framework of this model. The mixing 
parameter xd has been used to bound the intensity of these FCNC, still 
SM long distance contributions can be dominant [80 1. Automatically we 
get the short distance leading contribution to Br(D^ 1^^ |81|. The 
results are presented in flgure [20} as we can see, the natural range of the 
branching ratio for Br(D^ is of the order of 10~^^. Recently other 
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(a) Br(t ^ cZ) x 10^ (b) Br(t ^ iiZ) x 10^ 

Figure 19: Br(t cZ) and Br(t ^ uZ)^ Ax^ profiles in the sUVL model. 



surprising results in the charm sector have arised, as is the case of the direct 



CP violation in the channels tt+tt" and/or K^K' |82[]85] . The 

new tree level FCNC contribution has a suppression similar to the one we 
observe in Bt{D^ ji^ so we do not expect that sUVL FCNC can give 
an important contribution to this direct CP violation. Other mechanisms 
where vector-like quarks can play a role are at present under scrutiny |86|. 




Figure 20: Br(i^° ) and x^, Ax^ profiles in the sUVL model. 



7.3 Additional correlations and tut 

A very important issue in these models is the presence of a new non-sequential 
up type quark. This new quark decouples from the low energy theory as soon 
as its mass grows. This decoupling, in particular, refiects the fact that in this 



19 



limit all its couplings to light quarks vanish simultaneously. The consistency 
of the fit with this picture can be seen in the plots in figure [21] 




(c) rriT (GeV) vs. \VTb\- 

Figure 21: Correlations between tut and different \VTq\ in the sUVL model, 
68%, 95% and 99% CL regions (darker to lighter). 



We remind that recently DO [87|, CMS and ATLAS have published bounds 
on vector-like quark masses with values that range from 475 GeV in the CMS 
to 900 GeV in the ATLAS one (891. It has to be stressed that in 



case 



all these cases strong assumptions are made about the couplings of gauge 
bosons to this heavy quark and a light one. In particular CMS assumes that 
T ^ tZ 18 the dominant decay channel (100% branching ratio). In our case 
it is more likely to have a branching ratio around 25%, which reduces the 
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cross section by a substantial amount. In the ATLAS case it is implicitly 
assumed that \Vj^^\ ^ > 10"-^, that is at least a factor of ten larger than 
what flavour data allows; the corresponding cross section has to be reduced 
accordingly. It is out of the scope of this paper to analyse in detail the ttit 
lower bounds from direct production, but following the previous comments, 
it is clear that T masses as light as 350 GeV have to be considered. The 
most important result of our analysis is the likelihood profile for m^, shown 



in figure 22 
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Figure 22: ttit (GeV) Ax^ profile in the sUVL model. 



We can conclude that at 68% CL, tut < 0.6 TeV and at 95% CL tut < 
1 TeV. The reason of having upper bounds for these mass values is that 
vector like quarks decouple, so in the infinite singlet mass limit this model 
should reproduce the SM and cannot solve the known tensions. With a 
heavy new top T relatively light, in addition to an expected early production 
at the LHC, there will be other interesting signals: some of them will be 
also measured at the LHC. Several examples can be seen in the following 
plots. The oblique parameter AT controls figure [23| that clearly says that 
an important departure of from 1 will point to the presence of a light 
vector-like quark. The same happens with a sizeable departure from the SM 



value of several observables that we present in plots p4l p5^ and 26 



Aj/\^^^ A^j^ and A^^j^ — A^j^ can get values out of the SM range more 
easily for a relatively light T quark. Br(S+ t^i^t) and Br(S -^sT) can 
adjust better to its experimental values with a light m^. 
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Figure 23: ttit (GeV) vs. \Vtb\ in the sUVL model, 68%, 95% and 99% CL 
regions (darker to lighter). 
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Figure 24: tut (GeV) vs. Aj/^^ in the sUVL model, 68%, 95% and 99% CL 
regions (darker to lighter). 
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(a) rriT (GeV) vs. Br (5 X,7). 
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(b) rriT (GeV) vs. Br(5+ r+i/^). 



Figure 26: ttit (GeV) correlations in the sUVL model, 68%, 95% and 99% 
CL regions (darker to lighter) 
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8 Conclusions 



We have analysed most of the flavour data in the quark sector in the frame- 
work of one of the simplest extensions of the SM. We have incorporated an 
additional up quark of vector like character and singlet under SU{2)l. Its 
existence leads to the presence of FCNC and deviations of 3 x 3 unitarity, 
incorporated in a great variety of much more sophisticated models such as 
extra dimensions, little Higgs, etc 90-99). We have started by looking again 



to a model independent analysis to stress the fact that in the down sector 
the most important tensions appear in the bd sector. We realize that in the 
sUVL model, despite the violation of 3 x 3 unitarity, the situation is similar. 
In fact it is with the change of a unitary bd triangle to a quadrangle that 
we can improve the tensions among Br(5+ t^i^t) and ^j/^x^- By the 
same token, an enlargement of IV^^I - improving also Aj/^k^ - implies an 
enhancement of in the right direction of the DQ result of A^^^^ but not 
enough. Also ^^{Bd 1^'^ can get larger values compared to the SM 
when Br(5+ t^t^t) is in the right place. 

In the bs sector, in spite of the important constraints on mixing and CP 
violation coming from LHCb, it is still possible to have larger than SM values 
by some factors for A^^^ and Aj/\^^. Interestingly enough, larger values of 
Aj/\^^ are correlated with lower than SM values for Bt{Bs - but not 

too small -. In this respect, for this model, the rare decay BT{Bd 
seems more critical than the companion Bt{Bs ii^ It is also interesting 
to remark that if a deviation of the SM value of Aj/y^^ is found, then 7 
must be relatively large compared to its SM value. For the expected LHCb 
measurement, A^^^ — ^f^, we get natural values around (2 — 3) x 10~^, a 
factor of flve larger than the SM one. Needless to say, although we get a 
huge enhancement for still we are at more than 2a from the DO result. 
In the kaon sector, still an important enhancement of a factor of seven is 
possible in Bi^Kl n^vV). A very distintictive signature of this model is 
the enormous enhancement of the rare top decays t Zq (FCNC). It is 
t ^ Zc that gets the largest values, of the order of 10~^. Similar to these 
processes is Bi{D^ that can reach values of 5 x lO"-^-^. Finally 

we have to stress that the flt prefers a heavy top below 600(1000) GeV at 
68%(95%) CL. With this light up vector like quark we expect important 
deviations from the SM values of different observables. 

The flavour sector of the SM is being tested at an unprecedented level of 
accuracy and even more stringent tests will be available in the near future. 
In this paper we have studied how some of the present tensions between the 
SM and experimental data can be alleviated in the framework of a simple 
extension of the SM, where a up- type vector like quark is added to the SM. 
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Various correlations among physical observables are derived and provide a 
crucial set of experimental tests of the model. 
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A Experimental information 

This appendix summarizes the experimental information used in the analyses 
(see (Too -103]). For simplicity we present them in groups that share some 



important characteristic: observables related to (dominant) tree level decays, 
observables involving the mixing in different meson systems like B^s^ S^'s 
and kaons, rare decays of mesons and precision elect roweak information. 

Tree level observables are those whose extraction from experiment is pre- 
sumably unaffected by NP effects. These observables include moduli of the 
CKM elements in the first and second rows. Moduli of third row elements 
|V^^|, g = rf, 5, 6, are harder to extract. In fact the only relevant measurement 
is the one of the ratio of branching fractions R = Br(t Wb)/BT{t Wq)^ 
^ = + l^fsP + Notice that, when the mixing matrix 

is 3 X 3 unitary, this measurement is rather irrelevant since this unitar- 
ity requirement, together with the actual values of and iV^c^L forces 

|V^^| 1 — O{10~^). The picture changes when 3x3 unitarity is not as- 
sumed, since significant deviations (much larger than O{10~^)) of |V^^| from 
1 are possible. Finally, the physical phase 7 (in equation (T 
tree level observable. The actual values are collected in table 
T^Ur is also a tree level proces^ 

The next set of observables involves the mixings of B^ or Bg mesons (ex- 
cept for ^5'^, which involves both). We consider time-dependent CP asym- 



)), is also a 
The decay 



^Notice however that, as it is helicity suppressed and proportional to I^^^^P, NP sizeable 
contributions may appear in different beyond SM scenarios, but not in our case. 
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\vj 


0.97425 ± 0.00022 


\V 1 

\ us\ 


0.2252 ± 0.0009 


\vj 


0.230 ±0.011 


\v 1 

1 cs 1 


1.023 ±0.036 


\vj 


0.00389 ± 0.00044 




0.0406 ± 0.0013 


7 


(77± 14)° 


R 


0.88 ±0.07 


Br(5+ T+Ur) 


(1.68 ±0.31) X 10-4 





Table 1: Tree level observables [6H9l[T9l[20l[T0l] . 



metries Ajj^^s Ajj^^ (the "golden" channel in each system), mass and 
width differences AMb^^ AF^^, and AM^^, AF^, additional CP asymmetries 
involving different combinations of invariant phases, sin(2Q:), sin(2;5 + 7) and 
cos(2;5) (which removes a discrete ambigity in fixing 2^ = 8m~^{Aj/\iJKs))^ 
and, finally, semileptonic asymmetries A^^j^ and A^^j^. The actual val- 
ues are collected in table [2| We will also pay attention to the combination 
^SL ~ ^SL^ which LHCb results are expected. 





0.68 ± 0.02 




(0.508 ± 0.004) ps-i 




0.002 ± 0.0873 




(17.725 ± 0.049) ps'^ 


sin(2Q') 


0.00 ±0.15 


sm(2/3 + 7) 


1.0 ±0.16 


cos(2^) 


1.35 ± 0.34 




Ad 


-0.003 ± 0.0078 




-0.017 ±0.021 




-0.0017 ±0.0091 


AT, 


(0.116 ±0.019) ps-i 




-0.00787 ±0.00196 





Table 2: Bd and mixing-related observables (4[[5l[26l[27[ [TOOl[ToT] . 



Besides the observables involving the mixing, representative rare or sup- 
pressed decays are considered. The values in table [3] include the recent LHCb 
results used for Bt{Bs and Br(S^ that are playing a key 

role as the experimental bounds start to probe the SM expected range in the 
Bg case. The potential for discovery is still quite large in BT{Bd 
The actual values are collected in table |3l 

Moving from B mesons to kaons, representative observables that have 
to be considered address CP violation in K^-K^ mixing through ck and 
e' /ck-) and rare decays tt+z/P, Kl n^vV and Kl It is 

worth pointing that the experimental bound on n^vV is overseeded 

by a model independent bound, B>i{Kl ti^vV) < a Bt{K^ tt^uu) with 
a 4: A (see [79j). ji^ iJi~ is also expected to play a significant role since 

the absorptive part of the rate receives a contribution from an intermediate 
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Br{B ^ X,7) 


(3.56 ±0.25) X 10"^ 


Br(E X,ii+ii-) 


(1.60 ±0.51) X 10-6 


Br(5, ^ ^+^-) 


(0+25) X 10-9 







Table 3: Bd and 5, rare decays (7T}j75||TOO||ToT| . 



77 state that almost saturates the rate. The actual values are collected in 
table B 





(2.228 ±0.011) X 10-3 


s'/eK 


(1.67 ±0.16) X 10-3 




(1.73 ± 1.05) X 10-10 




< c(io-8) 


Bv{Kl ^ /i+A^-) 


(6.84 ±0.11) X 10-*^ 



Table 4: Kaon mixing and rare decays 1 78, 104-108 



For charm-related observables, we have D^-D^ mixing, for which we con- 
sider [Too[[ToT] 

xn = (0.8 X 10-2 . 

It is used to bound the short-distance contribution mediated by flavour 
changing Z-couplings. Notice that the experimental measurement is = 
(0.8 ±0.1) X 10~^. As in general D^-D^ mixing may receive substancial long- 
distance contributions, but the picture is quantitatively unclear, we content 
ourselves by requiring that the short-distance contribution does not overpro- 
duce this mixing (because in such a case one is implicitly resorting to long- 
distance contributions and a substantial cancellation to keep the theoretical 
expectation in place). Attention is also paid to the rare decay 
as GIM suppression yields a branching ratio below 10~^^ in the SM, while 
current bounds are at the 10"'' level an rates O{10~^^) are obtainable here. 

To complete the "menu" of relevant constraints, elect roweak precision 
data has to be considered. We do so through the oblique parameters T, 
and U 



109 



The T parameter is the most relevant one, while S is marginal 
and U negligible in the analysis. Their experimental values (with respect to 
the SM reference values), are [110] 



AS = 0.02 ± 0.11 , AT = 0.05 ± 0.12 
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with a correlation coefficient 0.879. 

The experimental values of the different observables quoted along this section 
are considered to follow gaussian profiles when the uncertainty is quoted in 
symmetric form, /xx±crx- This is extended to bounds by modelling them with 
gaussian profiles too, with central values equal to zero. The only observable 
not following that prescription is x^, for which no "sigma" is quoted for 
values below 0.8 x 10~^. As explained above, this only refiects the fact that 
we model values < 0.8 x 10~^ to be equally acceptable while a probability 
toll is paid for values > 0.8 x 10~^. 

Let us comment that, on the theoretical side, several predictions have a 
limited accuracy that cannot be neglected; this is the case, for example, of 
some QCD corrections and of several results produced through lattice compu- 



tations [111 -115 . They have associated uncertainties which are fundamental 
in the analyses. In particular some of the most important ones are quoted in 
table m 





1.243 ±0.028 




(275 ± 13) MeV 




(238.8 ± 9.5) MeV 




1.209 ±0.024 


Bk 


0.72 ±0.04 


Id 


(205 ± 9) MeV 



Table 5: Most relevant hadronic input. 



B Numerics 

In this appendix we explain the numerical procedure underlying the analyses 
and briefiy comment on statistics and the interpretation of the plots presented 
in the different sections. 

Considering the large amount of observables that are of potential interest 
in, (1), constraining the model or, (2), exploring the possibility of having pre- 
dictions experimentally distinguishable from SM prospects, an analytic ap- 
proach to use this information may be interesting for some particular scope, 
but falls short of the mark when the constraining power of different observ- 
ables (1) is not clear, (2) changes in different regions of parameter space or 
(3) when the number of available parameters is large. The strategy to deal 
with these difficulties relies on two cornerstones: considering enough observ- 
ables to overconstrain the new mixing matrix elements, and using numerical 
methods that allow an efficient exploration of the much larger parameter 
space that is available and drawing conclusions from it. The first ingredient 
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is at the heart of the selection of the set of observables hsted in appendix |A} 
Let us elaborate on the second one. For any given point in parameter space, 
we can compute the values of all the considered observables. Comparison of 
this set of predictions |44|[45l[50[|116 -118l with the corresponding measure- 
ments allows us to define, as usual, a likelihood function C (or equivalently a 

function = — 21n>C). This function has all the information we could be 
interested in: how good is the model to reproduce the experimental values of 
all the considered observables all over parameter space. The difficulty is that 
this function lives in a parameter space with too many dimensions to be easily 
grasped. This poses two problems: how do we explore this high-dimensional 
parameter space and how do we reduce the information to something that 
we can handle. The first one is a matter of the numerical methods available 
while the second one is a question of statistics. 

Numerical procedure. When dealing with such a number of observables 
and parameters, simple MonteCarlos are hardly appropriate for the task. 
The method of choice are Markov-chain MonteCarlos, in particular improved 
versions of the "classical" Metropolis algorithm. They allow an efficient ex- 
ploration of parameter space: as a tool for sampling the multidimensional 
likelihood function, obtaining probability distribution functions for a num- 
ber of interesting quantities is almost straightforward; in addition, through 
appropriate modifications, likelihood profiles are equally obtainable. 

Statistics . Reducing the information of jC in the n-dimensional param- 
eter space to zero, one or two-dimensional functions is a non-trivial task. 
There are two main "schools", bayesians and frequent ists. Without entering 
the polemic arena of the respective advantages and drawbacks we present 
frequent ist results: Ax^ profiles for individual quantities (from which confi- 
dence intervals can be easily obtained) and canonical 68%, 95% and 99% CL 
regions in two-dimensional plots. 

Beside the statistical treatment adopted, one may wonder if the overall 
goodness of a fit with the sUVL model compares well with the SM one: for 
that to be, the increase in the number of parameters of the model should 
improve significantly on the SM result. The quantity of interest for that 
comparison is the value of per degree of freedom for each fit. While for 
the SM we have 1.53 for 27 degrees of freedom, for the sUVL model we have 
1.56 for 21 degrees of freedom; the goodness of both fits is similar, i.e. the 
increase in the number of parameters in the sUVL really yields better fits to 
the existing data. 

One last comment closes this appendix. As far as the formulas we have 
used for the sUVL model, we have to stress that we have included for the 
first time in this type of analysis the right decoupling behaviour of the singlet 



up vector-like quark T 119-123 . In fact it amounts to use always the exact 
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formula for V — A coupling of the Z to up fermions {VV^)ij = 5ij — U-^Uj^ 
and never use the approximation {VV^)ij = Sij^ even if from a numerical 
point of view it may look irrelevant in some cases. 
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